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Abstract 

We formalize a classification of pair interactions based on the convergence properties of the forces 
acting on particles as a function of system size. We do so by considering the behavior of the proba- 
bility distribution function (PDF) P(F) of the force field F in a particle distribution in the limit that 
the size of the system is taken to infinity at constant particle density, i.e., in the "usual" thermody- 
namic limit. For a pair interaction potential V{r) with V(r — > oo) ^ l/r"' defining a bounded pair 
force, we show that P(F) converges continuously to a well-defined and rapidly decreasing PDF if 
and only if the pair force is absolutely integrable, i.e., for 7 > d— 1, where d is the spatial dimension. 
We refer to this case as dynamically short-range, because the dominant contribution to the force 
on a typical particle in this limit arises from particles in a finite neighborhood around it. For the 
dynamically long-range case, i.e., 7 < d — f , on the other hand, the dominant contribution to the 
force comes from the mean field due to the bulk, which becomes undefined in this limit. We discuss 
also how, for 7 < d — 1 (and notably, for the case of gravity, 7 = d — 2) P(F) may, in some cases, 
be defined in a weaker sense. This involves a regularization of the force summation which is gener- 
alization of the procedure employed to define gravitational forces in an infinite static homogeneous 
universe. We explain that the relevant classification in this context is, however, that which divides 
pair forces with 7 > d — 2 (or 7 < d — 2) , for which the PDF of the difference in forces is defined (or 
not defined) in the infinite system limit, without any regularization. In the former case dynamics 
can, as for the (marginal) case of gravity, be defined consistently in an infinite uniform system. 



PACS numbers: 98.80.-k, 05.70.-a, 02.50.-r, 05.40.-a 

I. INTRODUCTION 

Interactions are traditionally classified as long-range 
(or short-range) with respect to the non-additivity (or 
additivity) of the potential energy in the usual thermo- 
dynamic limit, i.e., when the number of particles N 
and volume V are taken to infinity at constant parti- 
cle density. This is the property which determines the 
way in which standard instruments of statistical mechan- 
ics are applied to determine equilibrium properties (see 
e.g. Indeed in the case of long-range interac- 

tions, these instruments are applied using an appropri- 
ately generalized thermodynamic limit, in which the cou- 
pling or density are also scaled with system size. Such 
an analysis gives rise generically to features at equilib- 
rium which are qualitatively different from those in short- 
range systems — inhomogeneous statistical equilibria, 
non-equivalence of statistical ensembles, negative specific 
heat in the microcanonical ensemble (see e.g. 0, Hi)- 
Most of these unusual features were first noted and stud- 
ied in the context of the study of gravitating systems 
in astrophysics (see e.g. [1, [1] for reviews), and it has 
been realized in recent years that they are more generic 
in long-range interacting systems. This thermodynamic 
analysis extended to long-range systems is believed to 



determine, however, the behavior of such systems only 
on time scales which diverge as some power of N (when 
expressed in terms of the characteristic dynamical time 
scales). On shorter times scales — usually those of in- 
terest in practical applications — study of several such 
systems (see e.g. |a-[lO| and references therein) shows 
that they appear generically, like in the well-documented 
case of gravity, to relax from almost any initial conditions 
to (almost) time-independent states — referred to vari- 
ously as "meta-equilibria" , "quasi-equilibria" or "quasi- 
stationary states" (QSS). The physics of these states, 
which are generically very different from those at ther- 
mal equilibrium, is understood to be the result of evolu- 
tion in the coUisionless regime described by Vlasov equa- 
tion (usually referred to as the "coUisionless Boltzmann 
equation" in the astrophysical literature [HI). Both the 
genesis of these states and their long-time relaxation are 
poorly understood, and are the subject of active study 
(see e.g. [lM3). 

In this article we consider a simple classification of 
pair interactions different to this usual thermodynamic 
one. Instead of considering the convergence properties 
of potential energy in the usual thermodynamic limit, we 
consider those of the force, in the same limit. The re- 
sulting classification can, like the usual one, be under- 



2 



stood easily from simple considerations. To see this let 
us consider, as illustrated schematically in Fig. 1, a uni- 
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FIG. 1: An (approximately) uniform system consisting of two 
sub-systems A and B. 

form system of particles interacting by a pair potential 
V{r — > oo) ^ l/r'^ , and divided into two pieces, A and B. 
For the usual thermodynamic classification one can con- 
sider that when the potential V{r) is integrable at large 
r, i.e., 7 > d where d is the spatial dimension, the poten- 
tial energy of a typical particle comes essentially from its 
interaction with particles in a finite region about it. The 
energy of a particle (e.g. P or P') in A is thus insensitive 
to whether B is present or not (and thus the total energy 
is equal, up to surface effects, to the sum of the energies 
of the subsystems): from an energetic point of view a 
particle "does not care" what the size of the system is, 
and the interaction is in this sense short-range. The dis- 
tinction we consider is the analogous one deduced when 
one reasons in terms oi force (or acceleration) rather than 
potential energy, and since forces are the primary physi- 
cal quantities in dynamics, we refer to the corresponding 
classification as one of dynamical range. It is straight- 
forward to see that in such a system, if the pair force 
is absolutely integrable, i.e., 7 > d — I, the force acting 
on a typical particle is due essentially to its interaction 
with particles in a finite neighborhood around it, while if 
7 < d — 1 this is not the case. Thus in the former case a 
particle in A "does not care" whether the sub-system B 
is present or not, and in this sense the interaction is "dy- 
namically short-range" . The classification differs from 
the standard one for interactions with d < 7 < d — 1 : for 
such interactions the potential energy "sees the bulk", 
but the force, which is its derivative, does not. 

While the principle motivation for defining such a clas- 
sification is that it may be relevant to understanding the 
qualitative behaviors of the out of equilibrium dynam- 
ics of such systems, it is not the aim of this article to 
establish that this is the case. We will limit ourselves 
in this respect to some brief remarks in our conclusions 
below. Our goal here is to provide a precise formula- 
tion of such a classification of the range of interactions 
based on the convergence properties of forces. While in 
the usual thermodynamic classification case one consid- 
ers (see e.g.[l[) the mathematical properties of essential 
functions describing systems at equilibrium in the limit 
iV — 00, F — 00 at fixed particle density Uq = N/V 



(i.e. the usual thermodynamic limit), we will consider 
the behavior of functions characterising the forces in this 
same limit. More specifically we consider, following an 
approach introduced by Chandrasekhar for the case of 
gravity [13, [ll] , the dcfinedness of the probability distri- 
bution function (PDF) of the force field in statistically 
homogeneous particle distributions as the size of the sys- 
tem becomes arbitrarily large. Such distributions and 
this limit may be described mathematically using the 
language of stochastic point processes, considering the 
class of such processes which have a well defined positive 
mean density when the infinite size limit is considered. 
To avoid any confusion we will refer to the usual thermo- 
dynamic limit in this context simply as the infinite sys- 
tem limit. Indeed the existence or non-existence of the 
quantities we are studying in this limit has no direct re- 
lation here to the determination of properties at thermal 
equilibrium. Further, in the context of the literature on 
long-range interactions the term "thermodynamic limit" 
is now widely associated with the generalized such limit, 
which involves adopting a different scaling of V (or possi- 
bly coupling constants) with N (for a discussion see e.g. 

In this article we also discuss, in Sec. lIVl a further (and 
different) classification which can be given of pair interac- 
tions according to their range. This is relevant when one 
addresses more generally, for any given pair interaction, 
a question which arises for Newtonian gravity in a cos- 
mological setting: can a consistent dynamics be defined 
in an infinite system with non-zero density? A rigorous 
approach to the same question and the connection with 
the possibility of defining a statistical mechanical state 
for the system has been developed in [2l| in the par- 
ticular case of a short range non-negative pair potential 
with finite support. Our conclusion, which generalizes a 
previous discussion given by two of us in [22| , is that the 
answer to this question is that a necessary and sufficient 
condition for such a dynamics is not the integrability of 
the pair force, but instead the integrability of its gradi- 
ent. This means that one requires 7 > d — 2, with gravity 
in any dimension (i.e. 7 = d — 2, the interaction potential 
solving the appropriate d-dimensional Poisson equation) 
being the marginal case in which such an infinite system 
limit may be defined. The reason is simply that, in an in- 
finite system without any preferred point (i.e. when this 
limit is defined respecting statistical translational invari- 
ance), the physically meaningful quantity is the relative 
position of particles (as there is no meaning to absolute 
position) . It is thus the convergence of relative forces on 
particles with system size which matters. In terms of the 
schema given above the distinction arises thus when one 
considers two close-by points (e.g. P and P' in Fig. 1) in 
sub-system A, say, and asks whether their relative forces 
— and thus relative motions — depend on the presence 
of B or not (or, equivalently, on the size of the system). 
The answer is that this difference of forces does not es- 
sentially depend on B if the gradient of the pair force is 
absolutely integrable, i.e., 7 > d — 2, as in this case this 
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difference is dominated by the contribution from parti- 
cles in a finite neighborhood around them. Thus for the 
case that d — 2<j<d— 1 the forces acting on two 
such particles become ill defined as the size of the sys- 
tem is extended to infinity, but their difference remains 
finite. Indeed, as has been discussed in [l^l in the context 
of gravity in one dimension, the diverging component of 
the force on a particle represents a force on their centre 
of mass, which has no physical relevance in an infinite 
system without a preferred origin. 



II. THE FORCE PDF IN UNIFORM 
STOCHASTIC POINT PROCESSES: GENERAL 
RESULTS 

We first recall the definitions of some basic quanti- 
ties used in the statistical characterization of a stochastic 
point process and define the total force PDF (see e.g. 
for a detailed discussion). We then derive some results 
on the analyticity properties of the latter quantity which 
we will exploit in deriving our central results in the next 
section. 



The paper is organized as follows. In the next section 
we recall the essential properties of stochastic point pro- 
cesses of relevance to our considerations, and then con- 
sider the general analyticity properties of the PDF of the 
total force at an arbitrary spatial point in such a parti- 
cle distribution. We show that, for any pair force which 
is bounded, this PDF in the infinite volume limit is ei- 
ther well defined and rapidly decreasing, or else vanishes 
pointwise, i.e., the total force is an ill defined stochas- 
tic quantity. This means that it suffices, when studying 
pair potentials with different possible behaviors at large 
scales, to show that some chosen moment of the PDF 
converges to a finite value in this limit (or diverges) in or- 
der to establish that the whole PDF itself is well-defined 
(or ill defined). In Sec. IIIII we give a general and for- 
mal expression for the variance of the total force PDF 
in a generic infinite uniform stochastic process in terms 
of the pair force and the two-point correlation proper- 
ties of the SPP. From this we then deduce our principal 
result that the force PDF exists strictly in the infinite 
system limit if and only if the pair force is absolutely 
integrable at large separations (i.e. 7 < d — 1) , while 
it can be defined only in a weaker sense, introducing a 
regularization, when the pair force is not absolutely in- 
tegrable. In the following section we discuss the physical 
relevance of the use of such a regularization, which is the 
generalization of a simple formulation given by Kiessling 
[23i] of that originally introduced by Jeans for the case of 
gravity 1241 . often misleadingly referred to as the "Jeans 
swindle" [ill [23I ] . By analyzing the evolution of density 
perturbations in an infinite system, we show that the 
physical relevance of such a regularization of the forces 
requires also a constraint on the behavior of the PDF of 
total force differences as a function of system size. This 
leads to the conclusion that 7 < d — 2 is the necessary 
and sufficient condition in order for it to be possible to 
have a well defined infinite system limit at constant den- 
sity for dynamics under a given pair interaction. In the 
conclusions we review briefly the the relation of our re- 
sults to previous work in the literature, and comment 
a little more on the possible relevance of our principle 
classification of interactions into dynamically short-range 
and dynamically long-range to the study of the out of 
equilibrium dynamics of such systems. 



A. Stochastic point processes 

In order to study the properties of the force field in 
the infinite system limit given by A'^ — > 00, — > 00 with 
fixed average density no > for a large scale uniform 
and spatially homogeneous particle system, we generalize 
the approach introduced by Chandrasekhar in for 
the total gravitational field in a homogeneous Poisson 
particle distribution to more general cases and spatial 
dimensions. To do so we need to characterize statistically 
point-particle distributions in this limit, and we do this 
using the language of stochastic point processes (SPP). 
The microscopic number density of a single realization of 
the process is 

n(x)=^J(x-xO (1) 

i 

where 6 is the d-dimensional Dirac delta function, x^ 
is the position of the i*'' system particle and the sum 
runs over all the particles of the system. We will limit 
our discussion to particle distributions in a euclidean 
d— dimensional space which are (i) statistically transla- 
tionally invariant (i.e. spatially homogeneous or station- 
ary) and (ii) large scale uniform in the infinite volume 
limit. Property (i) means that the statistical properties 
around a given spatial point of the particle distribution 
do not depend on the location of the point. In other 
words the statistical weights of two realizations of the 
point process, of which one is the rigidly translated ver- 
sion of the other, are the same and do not depend on the 
translation vector. In particular this implies that the en- 
semble average (i.e. average over the realizations of the 
SPP) (n(x)) of the microscopic number density takes a 
constant value hq > independent of x. Moreover the 
two-point correlation function of the microscopic density 
(n(x)n(x')) depends only on the vector distance x — x'. 
Feature (ii) means that the average particle number fluc- 
tuation 6N{R) = {{N^{R)) - (A^(i?))^)i/2 in a sphere 
of radius R increases slower with R than the average 
number {N{R))q V{R) with R, where V{R) oc R'^ is the 
volume of the d— dimensional sphere. 

Let us start by considering a generic realization of the 
particle distribution in a flnite volume V and let the total 
number of particles of the given realization be N. The 
particle positions x^ are fully characterized statistically 
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by the joint probability densty function (PDF) V^dxi}) 
conditional to having N particles in the realization ({x^} 
indicates the set of positions of all system particles in 
the given realization) . As a simple, but paradigmatic ex- 
ample we can think of the homogeneous d— dimensional 
Poisson point process. In this case VN{{xi}) = V^^ 
simply and independently of the value of uq. Given a 
function X{{xi}) of the N particle positions in thevol- 
ume V its average, conditional to the value of N, can be 
written as 



N 



N 



.1=1 



PAr({xa)X({x,}), 



where the position of each particle is integrated in the 
volume V. In order to evaluate the unconditional average 
of the property X, for which all possible outcomes of the 
value N are considered, one would need the probability 
qn of having N particles in the volume V, which permits 
to write: 



N ■ 



(2) 



^=0 



in a strict analogy with the grand canonical ensemble 
average in equilibrium statistical mechanics. However, 
since we are restricting the discussion to large scale uni- 
form particle distributions, for which 5N{R)/ {N{R)) 
vanishes for asymptotically large R, we expect that the 
larger the volume V the narrower will be the peak around 
N = {N{V)) = noV in which the measure qn will be con- 
centrated (for simplicity we have indicated with V both 
the region and its size). Asymptotically we expect that 
only the term of index NqV will contribute to the sum 
in Eq. i.e., for sufficiently large V we can write: 

(X) ~ {X)noV ■ 

In other words we can consider that for sufficiently large 
V the conditional PDF Vnovii^i}) characterizes com- 
pletely the statistical properties of the particle distribu- 
tion in the finite volume V and use this to evaluate in 
the following subsection the statistical properties of the 
total force. This is exactly what has been done, for in- 
stance, by Chandrasekhar in [l^ to calculate the total 
gravitational force PDF in the Poissonian case. 

In Appendix |X] we recall some of the basic definitions 
and properties of the statistical characterizations of uni- 
form SPP. We will use below notably two essential prop- 
erties of 5'(k), the structure factor (SF), which follow 
from its definition: 



B. General expression for the force PDF 

Let us consider now that the particles in any realization 
of the SPP interact through a pair force f (x), i.e., f (x) is 
the force exerted by a particle on another one at vectorial 
separation x. Further we will assume that the pair force 
is 



central, i.e. 



f (x) = x/(a;) , 



(5) 



where x x/x, and 

• bounded, i.e., there exists /q < oo such that |f (x)| = 
f{x) < fo for ah X. 

These assumptions simplify our calculations consider- 
ably, but do not limit our aim which is to establish the 
relation solely between the statistical properties of the 
force field and the behavior of the pair interaction at large 
distances. Note that the second assumption means that, 
in cases such as the gravitational or the Coulomb in- 
teraction, the divergence at zero separation is assumed 
appropriately regularized. We will briefly describe in our 
conclusions below how our results could be generalized 
to include such singularities. 

Let us assume for the moment that the system volume 
V is finite. As shown above, if V is sufficiently large, one 
can consider that the number of particles in this volume is 
deterministically NqV. We will deal with the important 
problem of the infinite volume limit defined by A^, — ^ oo 
with N/V — s> riQ > in the next subsection, by studying 
directly the limit V oo with fixed NqV. The total 
force field F(x) at a point x, i.e., the force on a test 
particle placed at a point x, may thus be written 



N 



N 



F(x)^^f(x-x,) = 5]. 



-/(|x-x,|). (6) 



The force field F(x) may be considered as a stochastic 
variable with respect to the SPP. Choosing arbitrarily 
the origin as the point where the total force is evaluated, 
the PDF of this force is formally defined by^ 



Piv(F) 



N 



II d^ 



Pw({x,})<5 



lim k'^Sik) 

fc-i-O 



0, 



i.e, the SF is an integrable function of k at fc = 0, 
and 



lim 5(k) = 1 



(3) ^ We consider here the unconditional force PDF, i.e., the force is 

that at an arbitrary spatial point, rather than that on a point 
occupied by a particle which belongs to the particle distribution. 
It is the latter case, of the conditional force PDF, which is often 
considered in calculations of this kind (see e.g. [25l427ll ). The 
distinction is not important here as the constraints we derive, 
which depend on the large scale correlation properties of the 

(4) particle distribution, would be expected to be the same in both 
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where we have used, as assumed, that f(— Xj) = — f(xi 
Using the identity 



1 



{2n) 



this can be rewritten as 



N 



f(. 



(7) 



The integral over the spatial coordinates in the above 
equation defines the characteristic function of the total 
field F 



Piv(q) 



N 



.4=1 



VN{{^^}) , (8) 



so that 



The integral over spatial configurations in Eq. ^ can 
be conveniently rewritten as an integral over the possible 
values of the pair forces due to each of the i = 1, ...,N 
particles: 



N 



.4=1 



QnHU}), (9) 



Pn (q) is an analytic function of q. Again from FT the- 
ory one has therefore that Pat (F) is a rapidly decreasing 
function of F: 

lim F°'PNiF) = 0, Va > 0. 

Thus Pn{F) is a well-defined function of which all mo- 
ments finite, i.e., < (|F|") < -|-oo for any n > 0. 

Let us now consider what happens when we take the 
limit — ;> oo with NoV. On one hand the joint PDF 
QAr({fi}) remains non-negative and absolutely intcgrable 
at all increasing V. On the other hand the support of this 
function remains compact with a diameter unaffected by 
the values of V, but fixed only by /q. Therefore we expect 
that the FT theorem keeps its validity also in the infinite 
system limit resulting in an analytical 

P(q) = Jim Pjv(q). 

N/Vo 

Therefore we will have that 

P(F) = hm PAr(F) 

Nov 



satisfies 



where 



lim F°'P(F) = , Va > 0. 

F— i-oo 



There are then only two possibilities for the behavior of 
Pn{(i) in the infinite system limit: 



QNim 



I 


' N 

l[d^x. 




.1=1 



N 



VN{{^^})Y[^[^^-f{^^)] 



(10) 

is the joint PDF for the pair forces f^. Note that, since F 
is the sum of the variables {fi} its characteristic function 
PAr(q) can be given as 



Piv(q) = Qjv({q. = q}) 



(11) 



where Q7v({qi}) is the A^d— dimensional FT of the joint 
pair forces PDF QAr({fi}), i.e.. 



Sjv({q.}) = 



N 



QNim). (12) 



C. Analyticity properties of the force PDF 

Prom the fact that the pair force is bounded it follows 
that QAr({fi}) has a compact support, and, since it is 
absolutely integrable (by definition), FT theory (see e.g. 
[28j ) implies that its characteristic function QnH^i}) is 
an analytic function of the variables {qi}. Consequently 



1. It converges to an absolutely integrable function 
which is not identically zero everywhere, giving a 
P{F) which is normalizable and non-negative on 
its support. Further all the integer moments of |F| 
are positive and finite. 

2. It converges to zero everywhere, giving P(F) = 
0. More specifically Pn{F) with NqV converges 
point-wise to the null function: it becomes broader 
and broader with increasing N (and V), but with 
an amplitude which decreases correspondingly and 
eventually goes to zero in the limit. 

This latter case is analogous to the case of the sum of 
identically distributed uncorrelated random variables: if 
this sum is not normalized with the appropriate power 
of the number A^ of such variables, the PDF of the sum 
vanishes point-wise in a similar way in the limit N ^ oo. 

In summary it follows from these considerations of the 
analyticity properties of -PAr(q) at increasing V that the 
case of a well defined, but fat tailed P(F), can be ex- 
cluded: in the infinite system limit the force PDF, if 
defined, is expected to be a normalizable and rapidly de- 
creasing function. 
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III. LARGE DISTANCE BEHAVIOR OF PAIR 
INTERACTIONS AND THE FORCE PDF 

In this section we use the resuh derived in the previous 
section to infer the main result of this paper: the relation 
between the large scale behavior of the pair interaction 
and the force PDF in the infinite system limit. 

We thus consider, as above, a central and bounded pair 
force such that 



We can thus write 



9 



.7+1 



for a; — > oo . 



(13) 



or, equivalently, a pair interaction corresponding to a 
two-body potential V{x) ~ g/{'jx'>') at large x for 7 / 
(and from V{x) ~ ^gliix for 7 — 0). Since the pair force 
is bounded, we have 7 > — 1. 

Given the final result derived in the previous section, 
it follows that, to determine whether the force PDF ex- 
ists, it is sufficient to analyze a single even moment of this 
PDF: because the PDF, when it exists, is rapidly decreas- 
ing, any such moment is necessarily finite and non-zero 
in this case, and diverges instead when the PDF does not 
exist. We choose to analyze the behavior of the second 
moment, (-F^), which is equal to the variance of the PDF 
since the first moment (F) is zero (see below) . We choose 
this moment because, as we will now see, it can be ex- 
pressed solely in terms of the FT of f (x) and of the SF of 
the microscopic density of the particle distribution. From 
these expressions we can then infer easily our result. 



A. Variance of the force in infinite system limit 

The formal expression of the total force acting on a test 
particle (i.e. the force field) at x in the infinite system 
limit may be written 



F(x) 



d'^x'-. 



-/(Ix-x'l)n(x') 



(14) 



where the integral is over the infinite space and n(x), 
given in Eq. ([T]), is the density field in a realization of 
the general class of uniform SPP we have discussed with 
positive mean density uq. 

It is simple to show, using Eq. and the definition 
of the SF given above in Eq. (jAl[) , that formally 



1 



{2n) 



d''k\i{k)\^S{k) 



(15) 



where f(k) is the (d-dimensional) FT of xf{x). It is 
straightforward to show that f(k) = k/(fc), where the 
explicit expression for f{k) is given in the appendix^. 



^ Note that only in d = 1 does f{k) coincide with the direct FT of 



(16) 



dkk'^-^\f{k)YS{k) , 



2d-i7rd/2r(d/2) 7o 
where T{x) is the usual Euler Gamma function. 

B. Force PDF for an integrable pair force 

Let us now consider the integrability of the integrand 
in Eq. (fT6|) . We start with the case in which f{x) is not 
only bounded but integrable in E'', i.e., with 7 > d — 1. 
Given these properties, it is straightforward to verify, us- 
ing the conditions ([3]) and (HJ on S{k) and standard 
FT theorems, that the function |/(fc)p5(fc)) is also inte- 
grable in M''. The variance is therefore finite, from which 
it follows that the PDF exists, and furthermore that all 
its moments are finite. 

C. Force PDF for a non-integrable pair forces 

For a pair force which is absolutely non-integrable, i.e., 
7 < d - 1, the FT f(k) of f(x) in Eq. ((16]) is defined only 
in the sense of distributions, i.e., the integrals over all 
space of f{x) must be defined by a symmetric limiting 
procedure. Physically this means that the expression Eq. 

for the force on a particle in infinite space must be 
calculated as 



F(x) = hm hm 

/j->0+ V-5.0. 



-/(|x-x'|)e-^l''-'''l7i(x')dV, 



X X 

(17) 

where the two limits do not commute. In other words, 
F(x) is defined as the zero screening limit of a screened 
version of the simple power law interaction in an infinite 
system. The expression Eq. p6)) is then meaningful when 
/(fc) is taken to be defined in the analogous manner with 
the two limits /i 0+ of the screening and V ^ 00 (i.e. 
with the minimal non-zero mode k ^ l/V — > 0+) taken 
in the same order as indicated in Eq. (|17p . 

Let us consider then again, for the case 7 < d — 1, the 
integrability of the integrand in Eq. (fT6)) . To do so we 
need to examine in detail the small k behavior of /(fc). It 
is shown in the appendix that, as one would expect from a 
simple dimensional analysis, for /(r 00) we 
have f{k 0) k^'^^''^^ in any d, for the case of a pair 
force which is not absolutely integrable, and bounded, 
i.e., — l<7<d— 1. It follows then from Eq. (fT6)) that 
the variance is finite for a given 7 only for a sub-class of 
uniform point processes, specifically those which satisfy 



lim k 

k^O 



-rf+27+2 



S'(fc) = 0, 



i.e., for S{k ^ 0) - fc" with 

n > d - 27 - 2 = -d + 2(d - 1 - 7) . 



(18) 



(19) 
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For uniform point processes violating this condition, i.e., 
with S{k -J- 0) fc" and -d < n < -d + 2{d--f-l), the 
variance diverges. It follows from the results on the PDF 
of F presented in the previous section that the total force 
itself F(x) is then badly defined in the infinite system 
limit. 

These results of Sec. IIII Bl and Sec. IIII CI combined are 
the central ones in this paper, anticipated in the intro- 
duction. 

Firstly, when pair forces are absolutely integrable at 
large separations, the total force PDF is well defined in 
the infinite system limit, while for pair forces which are 
not absolutely integrable this quantity is ill defined. This 
has the simple physical meaning anticipated in the intro- 
duction: when this PDF is well defined, the force on 
a typical particle takes its dominant contribution from 
particles in a finite region around it; when instead the 
PDF is ill defined far-away contributions to the total 
force dominate, diverging with the size of the system. 
Thus absolutely integrable pair forces with 7 > d — 1 
are, in this precise sense, "short-range", while they are 
"long-range" when 7 < fi — 1. To avoid confusion with 
the usual classification of the range of interactions based 
on the integrability properties of the interaction poten- 
tial, we will adopt the nomenclature that interactions in 
the case 7 > d — 1 are dynamically short-range, while 
for 7 < d — 1 they are dynamically long-range. Thus an 
interaction with d — 1 < 7 < d can be described as ther- 
modynamically long-range but dynamically short-range. 

Secondly the results in Sec. IIII CI detail how, for 7 < 
d — 1, the force PDF in the infinite system limit may 
be defined provided an additional prescription is given 
for the calculation of the force. In the next section we 
explain the physical meaning and relevance of this result. 



IV. DEFINEDNESS OF DYNAMICS IN AN 
INFINITE UNIFORM SYSTEM 

The regularization Eq. (1171) is simply the generalization 
to a generic pair force with 7 < d — 1 of one which is 
used for the case of Newtonian gra vity , often referred to 
as the "Jeans swindle" (see e.g. [ll|). It was indeed 
originally introduced by Jeans [IJl in his treatment of 
self-gravitating matter in an infinite universe. However, 
as explained by Kiessling in [23| . its denomination as 
a "swindle" is very misleading, as it can be formulated 
in a mathematically rigorous and physically meaningful 
manner, precisely as in Eq. (|17p . 

The prescription Eq. (|17p simply makes the force on a 
particle defined by setting to zero the ill defined contri- 
bution due to the non-zero mean density: 

(F(x))- hm no / ii^/(|x-x'|)e-^l---'ldV - , 
(,^0+ J |x-x'| 

(20) 



The force on a particle can thus be written as 

F(x)= lim^y ^_^/(|x_x'|)e-'^l"-"'l<5n(x')dV, 

(21) 

where 6n{x') = n(x') — uq is the density fluctuation field. 
It is straightforward to show that the derived constraint 
(|19|) corresponds simply to that which can be antici- 
pated by a naive analysis of the convergence of the inte- 
gral Eq. (mj: treating (5n(x') as a deterministic function 
(rather than a stochastic fleld) one can require it to decay 
at large |x'| with a sufficiently large exponent in order to 
give integrability; taking the FT to infer the behavior of 
(5n(k)p one obtains the condition 

The relevance of the results we have derived for the 
force PDF in the inflnite system limit using this regular- 
ization arises thus, as it does in the case of Newtonian 
gravity, when one addresses the following question: is it 
possible to define consistently dynamics under a given 
pair interaction in an infinite system which is uniform at 
large scales? As we now discuss, generalizing consider- 
ations given by two of us in 22] for the specific case of 
gravity in d = 1, the answer to this question is in fact 
phrased in terms of the definedness of the PDF of force 
differences rather than that of forces. This leads then to 
our second classification of pair interactions. 



A. Evolution of fluctuations and definedness of 
PDF 

Let us consider first an infinite particle distribution 
which is such that the total force PDF is defined at some 
given time, i.e., for 7 > d — 1 we may consider any uni- 
form SSP, while for 7 < d— 1 we may consider (employing 
the regularization discussed) only the class of SSP with 
fluctuations at large scales obeying the condition ([TOl) at 
this time. The forces on particles at this initial time are 
then well defined. This will only remain true, however, 
after a finite time interval, if the evolved distribution con- 
tinues to obey the same condition . Let us determine 
when this is the case or not. 

In order to do so, it suffices to consider the evolution 
of the density fluctuations, and specifically of the SF at 
small k, due to the action of this force field. Given that 
we are interested in the long-wavelength modes of the 
density field, we can apply the differential form of the 
continuity equation for the mass (and thus number) den- 
sity between an initial time t = and a time t = St: 

n(x, dt) - n(x, 0) = V[n(x, 0)u(x, 0)] (22) 

where u(x, 0) is the infinitesimal displacement field. Sub- 
tracting the mean density no from both sides, and lin- 
earizing in SnipCjSt) = [n(x, (5i) — no] smd u(x, 0), we 
obtain, on taking the FT, 

Sn{k, 5t) = Sn{k, 0) -1- i no k • u(k, 0) . (23) 
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Taking the square modulus of both sides, in the same 
approximation we get 

\Sn{k,St)f - |^n(k,0)|2 = (24) 
n^fc2|u(k)p + 2knolm[5n(k,0)u*(k,0)] . 

If the displacements are generated solely by the forces 
acting (i.e. assuming velocities are initially zero), we have 
that 

u(x,0) = iF(x, 0)5*2 (25) 

and thus, that |u(k)p oc |F(k)|2. The latter quantity is 
given, using Eq. ([T5|). by 

|F(k)p = |/(fc)|2^(fc). (26) 

In the analysis in the previous section we used the result 
that at smah fc, /(fc) ~ fc-'^+T+i. Thus |u(k)|2 ~ fc2™+"^ 
where m = — (i + 7+ 1, if S'(fc) ^ k". It then follows, 
from Eq. (|24l) . that the small k behavior of the time- 
evolved SF is given by 

Sst{k ^ 0) ~ fc" + fcl+™+« + fc2+2m+n (27) 

It can be inferred that the leading small k behavior of the 
SF is unchanged if and only if m + 1 > 0, i.e., 7 > d — 2. 
Gravity (7 = d — 2) is the marginal case is which the long 
wavelength contribution to the SF generated by the evo- 
lution has the same exponent as the initial SF: this is the 
well known phenomenon of linear amplification of initial 
density perturbations (see e.g. f29] ) which applies'^ 
in infinite self-gravitating systems (derived originally by 
Jeans). 

If, on the other hand, 7 < d — 2 (i.e. the interac- 
tion is "more long-range" than gravity in d dimensions) 
the exponent of the small k behavior is reduced from 
n to n — 2((i — 2 — 7). Given that our result is for an 
infinitesimal time 5t, this indicates in fact a pathologi- 
cal behavior: in any finite time interval the exponent n 
should become, apparently, arbitrarily large and nega- 
tive, while, as shown in Sect. |TT1 the constraint n > —d 
is imposed by the assumed large scale uniformity of the 
SPP. In other words this result means that, in the infi- 
nite system limit, when 7 < d — 2, the condition of large 
scale uniformity is violated immediately by the dynami- 
cal evolution. The reason is simply that in this case the 
rate of growth of a perturbation at a given scale increases 
with the scale. Indeed this is the essential content of the 
analysis given just above: through the continuity equa- 
tion, the perturbation to the density field is proportional 



^ The result does not apply, however, when n > 4 j29|; the reason is 
that fluctuations with S{k — > 0) ~ fc'' arise generically from any 
rearrangement of matter due to dynamics which conserves mass 
and momentum locally. These effects are neglected implicitly 
above when we use the continuum approximation to the density 
fluctuation field. 



to the gradient of the displacement field, which in turn is 
simply proportional to the gradient of the force. As we 
now detail more explicitly , when 7 < d — 2, this quantity 
diverges with the size of the system. 



B. PDF of force differences 

Let us consider now the behavior of the PDF of the dif- 
ference of the forces between two spatial points separated 
by a fixed vector distance a: 

AF(x;x + a) = F(x) - F(x + a) . (28) 

If this quantity is well defined in the infinite system limit, 
its PDF 'P{A.F; a) will be independent of x and will have 
a parametric dependece only on a = |a| because of the as- 
sumed statistical translational and rotational invariance 
of the particle distribution. 

The analysis of the properties of ^(AF; a) in the in- 
finite volume limit is formally exactly the same as that 
given above for the total force F, with the only replace- 
ment of the pair force in Eq. (|13p by the pair force dif- 
ference: 

Af (x, X + a) = f (x) - f (x -I- a) , (29) 

i.e., the difference of the pair forces on two points lo- 
cated at X and x -t- a due to a point at the origin. As- 
suming again the possible small scale singularities in this 
pair force difference to be suitably regulated, our previ- 
ous analysis carries through, the only significant change 
being that, as a; -^■ cx), 

Af(x,x + a) ~ax/a;'^+^ (30) 

Proceeding in exactly the same manner to analyse 
V{AF;a), we find that 

• For 7 > d — 2, i.e., if the gradient of the pair force 
at fixed a is an absolutely integrable function of x 
at large separations, the PDF ^^(AF; a) is well de- 
fined in the infinite system limit, and is a rapidly 
decreasing function of its argument for any SPP. 
This is true without any regularization. 

• For 7 < d—2, on the other hand, a well defined PDF 
may be obtained only by using the regularization 
like that introduced above in Eq. (fTT]) . Therefore 
the PDF of the force differences then remains well 
defined, i.e., the force difference AF(x;a) remains 
finite at all x, only in a sub-class of SPP defined by 
the constraint 

n>d-27--4 = -d + 2(d-2-7). (31) 

For the case of gravity 7 = d — 2 this coincides with 
the full class of uniform SPP, while for any smaller 
7, it restricts to a sub-class of the latter. 
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C. Conditions for definedness of dynamics in an 
infinite system 

Our analysis in Sec. IIV Al of the evolution of density 
perturbations under the effect of the mutual pair forces 
gave the sufficient condition 7 > d — 2 for the consistency 
of the dynamics in the infinite system limit, but with the 
assumption that the total force PDF was itself defined. 
This means that, in the range d—2 < 7 < d—1, the result 
derived applies only to the sub-class of infinite uniform 
particle distributions in which the large scale fiuctuations 
obey the condition (fT9)) . It is straightforward to verify, 
however, that the analysis and conclusions of Sec. IIV Al 
can be generalized to cover all uniform SPP for 7 > d — 2. 
In line with the discussion given above, the analysis re- 
quires in fact only assumptions about the behavior of 
the gradient of the forces, rather the forces themselves. 
More specifically, the only equation which explicitly con- 
tains the force, Eq. (j25p . is a purely formal step which can 
be modified to include the possibility that the force di- 
verges with system size. Indeed if the force — at a given 
point — includes such a divergence it is sufficient that 
this divergence cancels out when we calculate the differ- 
ence between this force and that at a neighboring point. 
Physically this means simply that, as discussed above, 
when we consider the relative motions of particles, it is 
sufficient to consider relative forces. Further, as we are 
considering the limit of an infinite system in which there 
is no preferred point (i.e. statistical homogeneity holds), 
only relative motions of points has physical significance, 
and therefore only the spatial variation of the forces can 
have physical meaning. These latter statements can be 
viewed as a kind of corollary to Mach's principle: if the 
mass distribution of the universe is, as it is in the case we 
consider, such that there is no preferred point in space 
(and, specifically, no center of mass) inertial frames which 
give absolute meaning to forces (rather than tidal forces) 
cannot be defined. 

In summary our conclusion is that the necessary and 
sufficient condition for dynamics to be defined in the in- 
finite system limit — in analogy to how it is defined for 
Newtonian self-gravitating particles in a infinite universe 
of constant density — is that the gradient of the pair force 
be absolutely integrable at large separations. Gravity is 
the marginal (logarithmically divergent) case in which 
such a dynamics can be defined, but only by using a pre- 
scription such as Eq. (IT71) . Further these conditions on 
the range of pair forces can be expressed simply as one 
on the existence of the PDF of force differences of points 
as finite separations in the infinite system limit. 



V. DISCUSSION AND CONCLUSIONS 

In conclusion we make some brief remarks on how the 
results derived here relate to previous work in the lit- 
erature on force PDFs. In this context we also discuss 
the important assumption we made throughout the arti- 



cle, that the pair force considered was bounded. Finally 
we return briefly to the question of the relevance of the 
classification dividing interactions according to the inte- 
grability properties of the pair force, concerning which 
we have reported initial results elsewhere J3(| . 

The first and most known calculation of the force PDF 
is that of Chandrasekhar [13], who evaluated it for the 
gravitational pair interaction in an infinite homogeneous 
Poisson particle distribution (in d — 3). This results in 
the so-called Holtzmark distribution, a probability distri- 
bution belonging to the Levy class (i.e. power law tailed 
with a diverging second moment) with P(F) ~ p~^/^ 
at large F. According to our results here, a well de- 
fined PDF may be obtained for such a force law, which 
is not absolutely integrable at large separations, only 
by using a prescription for the calculation of the force 
in the infinite system limit. In his calculation Chan- 
drasekhar indeed obtains the force on a point by sum- 
ming the contributions from mass in spheres of radius 
R centered on the point considered, and then taking 
i? — >■ 00 (with no fixed). This prescription is a slight 
variant of the one we have employed (following Kiessling 
[2^): instead of the smooth exponential screening of the 
interaction, it uses a "spherical top-hat" screening so 
that the force may be written formally as in Eq. (jI7p 
with the replacement of e^^'''^'' ' by a Heaviside func- 
tion Q{fJ,~^ — |x — x'l). It is straightforward to verify that 
the result of Chandrasekhar is unchanged if the smooth 
prescription Eq. ([T7)) is used instead. As the Poisson 
distribution corresponds to an SF S{k — ?> 0) ~ fc" with 
n = 0, the general condition for the existence of the 
PDF we have derived, which gives n > —I for gravity 
in d = 3, is indeed satisfied. The fact that the PDF is 
power-law tailed (and thus not rapidly decreasing) arises 
from the fact that the calculation of Chandrasekhar does 
not, as done here, assume that the singularity in the grav- 
itational interaction is regularized. Indeed it is simple to 
show explicitly [T^ that this power law tail arises from 
the divergence in the pair force at zero separation. This 
can be done by considering the contribution to the to- 
tal force on a system particle due to its nearest neighbor 
particle, which turns out to have a power law tail identi- 
cal, both in exponent and amplitude, to that of the full 
P(F). 

Our analysis shows that it is true in general that well 
defined, but power-law tailed force PDFs, can arise only 
when there are singularities in the pair force: for a 
bounded force we have seen that the PDF is necessar- 
ily rapidly decreasing when it exists. More specifically, 
returning to the analysis of Sec. Ill C[ it is straightforward 
to see that the crucial property we used of QAr({fi}), that 
it have compact support, is no longer valid when the pair 
force has singularities. The analyticity properties which 
lead to a rapidly decreasing PDF may then not be in- 
ferred. We note that this is true at finite N, and has 
nothing to do with the infinite volume limit, i.e., the 
appearance of the associated power-law tail arises from 
the possibility of having a single particle which give an 
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unbounded contribution rather than from the combina- 
tion of the contribution of many particles which then 
diverges in the infinite system hmit. The exponent in 
such a power-law tail will depend on the nature of the 
divergence at small separation. More specifically, for a 
central pair force as considered above and now with a 
singularity f(x — > 0) ~ 1/x", a simple generalization 
of the analysis for the case of gravity (see 18]) of the 
leading contribution to the total force coming from the 
nearest neighbor particle leads to the conclusion that 
P{F ^ cx)) ~ F-'^-^ (where F = |F|). This implies 
that the variance diverges (i.e. the PDF becomes fat- 
tailed) for a > d/2. 

Force PDFs have been calculated in various other spe- 
cific cases. Wesenberg and Molmer [2^ derived that of 
forces exerted by randomly distributed dipoles in d = 3, 
corresponding to a pair force with j — 2. According to 
our results this is the marginal case in which a summation 
prescription is required for the force, and indeed a pre- 
scription using spheres, like that used by Chandrasekhar 
for gravity, is employed. We note that [25^ focusses on 
the power-law tails associated with the singularity at zero 
separation of the force, which lead in this case (as can be 
inferred from the result summarized above) to the diver- 
gence of the first moment of the force PDF. One of us 
(AG) has given results previously [1^ for the PDF for 
a generic power-law interaction in d = 1 for 7 > — 1 
in our notation above. The conditional force PDF is 
then derived for the case of an infinite "shuffled lattice" 
of particles, i.e., particles initially on an infinite lattice 
and then subjected to wncorreZated displacements of finite 
variance, and using again, as Chandrasekhar, a "spher- 
ical top-hat" prescription for the force summation (for 
7 < 0, when the pair force is not absolutely integrable). 
It is simple to show fisl that such a distribution has an 
SF with n = 2 at small k, and thus the existence of the 
force PDF in these cases is again in line with the con- 
straint derived. Power-law tails are again observed 
in these cases, and their exponents related explicitly to 
the singularity in the assumed power-law force at zero 
separation. 

The calculation of Chandrasekhar has been general- 
ized in [l^l to the case of particles on an infinite shuffled 
lattice. This leads again, in line with condition ([T^. 
to a well defined PDF, again with or without power-law 
tails according to whether the singularities in the pair 
force are included or not. Chavanis [3l| considers, on the 
other hand, the generalization of Chandrasekhar calcu- 
lation (for the PDF of gravitational forces in a Poisson 
distribution) to c? = 2 and d = 1. The condition for 
gravity (7 = c? — 2) gives n > —d + 2, which implies that 
the force PDF is not defined in the infinite system limit 
we have considered for d < 2, and indeed in 31] well de- 
fined PDFs are obtained in d = 2 and d = 1 by using a 
different limiting procedure involving in each case an ap- 
propriate rescaling of the coupling with N. T he p hysical 
meaning of such a procedure is discussed in |22| . which 
considers in detail the calculation of the force PDF for 



gravity in d = 1 in a Poisson distribution (as in 3l|). An 
exact calculation of the force PDF of the screened gravi- 
tational force in the infinite system limit is given, which 
allows one to see in this case exactly how the general 
result given here is verified in this specific case: all mo- 
ments of the PDF diverge simultaneously as the screening 
length is taken to infinity, giving a PDF which converges 
point-wise to zero. The force PDF for gravity in d = 1 for 
a class of infinite particle distributions generated by per- 
turbing a lattice has been derived recently by three of us 
in . It is straightforward to show that one of the con- 
ditions imposed on the perturbations to obtain the PDF, 
that the variance of the perturbations be finite, corre- 
sponds in fact to the condition n > 1 which coincides 
precisely with the more general condition (|19p derived 
here. Unlike in the other specific cases just discussed, it 
turns out that in this case (gravity in d = 1) it is in fact 
necessary to use the smooth prescription Eq. (fTT]) . As 
explained in detail in (32| . the top-hat prescription does 
not give a well defined result in this case, because sur- 
face contributions to the force which do not decay with 
distance in this case are not regulated by it. We under- 
line that the general result given in the present article 
are for this specific prescription Eq. (|17p. Further anal- 
ysis would be required to derive the general conditions 
in which a top-hat prescription also gives the same (and 
well-defined) PDF. 

Finally let us comment on why we anticipate the classi- 
fication of pair interactions according to their "dynamical 
range", formalized here using the force PDF, should be 
a useful and relevant one physically in the study of sys- 
tems with long-range interactions. The reason is that this 
classification refiects, as we have explained, the relative 
importance of the mean field contribution to the force 
on a particle, due to the bulk, compared with that due 
to nearby particles. Now it is precisely the domination 
by the former which is understood to give the regime of 
coUisionless dynamics which is expected to lead to the 
formation of QSS states, which are usually interpreted 
to be stationary states of the Vlasov equations describ- 
ing such a regime of the dynamics (see e.g. [HI). In a 
recent article [30| by three of us, we have reported a nu- 
merical and analytical study which provides strong ev- 
idence for the following result, very much in line with 
this naive expectation: systems of particles interacting 
by attractive power law pair interactions like those con- 
sidered here can always give rise to QSS; however when 
the pair force is dynamically short-range their existence 
requires the presence of a sufficiently large soft core, while 
in the dynamically long-range case QSS can occur inde- 
pendently of the core, whether hard or soft, provided it 
is sufficiently small. In other words only in the case of 
a pair force which is "dynamically long-range" can the 
occurrence of QSS be considered to be the result only of 
the long distance behavior of the interaction alone. This 
finding is very consistent with what could be anticipated 
from the preceding (naive) argument: the effect of a "soft 
core" is precisely to reduce the contribution to the force 
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due to nearby particles, which would otherwise dominate 
over the mean field force in the case of a pair force which 
is absolutely integrable at large distances. Indeed the 
meaning of "sufficiently large" specified in [s^ is that 
the size of the soft core must increase in an appropriate 
manner with the size of the system as the limit N ^ oo 
is taken, while we have always implicitly assumed it to be 
fixed in units of the interparticle distance here. Further 
work on these issues will be reported elsewhere. 

We thank M. Kiessling and T. Worrakitpoonpon for 
useful conversations, and P. Viot for useful comments on 
the manuscript. 



Appendix A: One and two point properties of 
uniform SPP 



In this appendix we give the general one and two-point 
statistical characterization of a SPP which is uniform on 
large scales. 

The description of the correlation properties of a 
generic uniform SPP is given by the 71-point correla- 
tion functions of the density field. For our considera- 
tions it will turn out to be sufficient to consider only 
the two-point properties, and more specifically it will be 
most convenient to characterize them in reciprocal space 
through the structure factor (SF) (or power spectrum). 
This is defined by 



where 



5(k) 



Sn(k; V) 



lim 

V— >CXD 



\Snik;V)\' 
noV 



■[n(x) - no] 



(Al) 



(A2) 



With these normalisations the SF of an uncorrelated 
Poisson process is S(k) = 1. For a statistically isotropic 
point process S'(k) = S{k), where k = |k|. We recall here 
that S'(k) is the Fourier transform (FT) of the connected 
two point density correlation function: 



5(k) = J d'^a 



-ik-x 



where 
C(x) 



(n(xo +x)n(xo)) - 



no 



C(x) 



= (5(x) + no/i(x) 



In the last expression we have explicitly separated in 
the correlation function C(x) the shot noise term i5(x), 
present in all SPP and due to the "granularity" of 
the particle distribution, from the "off-diagonal" term 
noh(yi) which gives the actual spatial correlations be- 
tween different particles. 

In the paper we study the convergence properties of 
forces at large distances and are thus mainly interested 



in the properties of the SF at small k. In this respect we 
will use the following limit on the SF which follows from 
the assumed uniformity of the SPP: 

lim k'^SCk) = , 

fe^O 

i.e, the SF is an integrable function of k at fc = 0. This 
constraint simply translates in reciprocal space the re- 
quirement from uniformity on the decay of relative fluc- 
tuations of the number of particles contained in a volume 
V about the mean at large V: 



lim 



{NjVr) - {N{V)y 
{N{V))^ 



0. 



Given that {N{V)) oc V, the root mean square fluctua- 
tion of particle number A'^ in a volume V must diverge 
slower than the volume V itself in order that this condi- 
tion be fulfilled. (This is equivalent to saying that C(x) 
must vanish at large x). 

We use likewise in the paper only one constraint on 
the large k behavior of the SF, which is valid for any 
uniform SPP (see e.g. [l^l) and coincides with the shot 
noise term in the correlation function C(x): 

lim S'(k) = 1 . 



Appendix B: Small k behavior of f(k) 

We are interested in the small k behavior of the Fourier 
transform f(k) of the pair force in d dimensions in the 
case where the pair force f(x) = x/(a;), where x = 
is non-integrable but converges to zero at a; — > 00, i.e., 
/(r) ^ at large x with — 1 < 7 < d — 1. 

We first show that for a function f(x) = xf{x), its 
Fourier transform, f(k) = FT[f(x)](k), can be written 
f(k) = k V'(fc) where ip{k) is a function depending only 

In order to obtain this 



on the modulus of k and k = t}^ 
result, we start by writing 



f(k) 



d'^x f(x) 



d'^x xf{x)e~ 



ikx 



where this integral is defined in the sense of functions or 
distributions according to the integrability of f{x). 
In the following we denote by (§1, §2, . . . , e„) the 
cartesian vector basis in d-dimension and we define 
(r, 01,02, ■■■ , Od-i) the hyper-spherical coordinates of x. 
Considering k = fc ei and denoting for simplicity 6 — 61, 
we can write 



where 



f(k) = J d'^x if{x) 



a— ± 

^'^^ = ( n sin^{Sd-j)ded-j^x'^^^dx . 
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Projecting f(k) on the cartesian basis, it is easy to see 
that the only non- vanishing term is ei.f(k) which gives 



the small k behavior is determined by this leading diver- 
gence at a; — i> oo, 



3i./(k) = Ce,_^i / dxx 



„d-l 



X d9 sin''- \e) cos ef{x) 6-'"''''°'"', 
Jq 

where Cg-^-^ is a constant term coming from the inte- 
gration over all the hyper-spherical coordinates 9i with 
i 1. We thus can write f(k) = k ip{k) where ilj{k) is a 
function depending only on the modulus of k. 

We now focus our attention on the small k behavior of 
the term 



dxx'^-\f{r)e 



-ikxcosO 



(Bl) 



where the function f{x) is non-integrable but converges 
to zero at a; — oo, i.e., /(x) x"^'*'^^' at large x 
with —1 < 7 < d — 1, and thus can be written f{x) — 
2;-(7+i) h{x) with h{x) a smooth function, integrable 
at X = and such that x^+^/i(x) — > for x — >■ oo. 
Defining explicitly ea. (|Bip in the sense of distributions. 



lim / dxx'^-^ ^-^kxcost 



x-^+i 



(B2) 



where the parameter /i > 0. We define a = c? — 7 — 2 
which satisfies —l<a<d—l and rewrite eq. (|B2[) 



lim / dxx"e-(*'=™"^+^)"= . 

This can be easily calculated with Laplace's transform 
and gives 



+ ifccos6')"+i 



We can conclude that 



lim / dxx'^-^— -re-*'^^™''^ 



r("+i) cos-("+i)(6l)r(a + l)A:^("+i) ~ k^'"^+^ 
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